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ABSTRACT 

We undertake an investigation into the spatial structure of the M31 satellite system utilizing the distance 
distributions presented in a previous publication. These distances make use of the unique combination of depth 
and spatial coverage of the Pan- Andromeda Archaeological Survey (PAndAS) to provide a large, homogeneous 
sample consisting of 27 of MSl's satellites, as well as M31 itself. We find that the satellite distribution, when 
viewed as a whole, is no more planar than one would expect from a random distribution of equal size. A disk 
consisting of 15 of the satellites is however found to be highly significant, and strikingly thin, with a root-mean- 
square thickness of just 12.34^^43 kpc. This disk is oriented approximately edge on with respect to the Milky 
Way and almost perpendicular to the Milky Way disk. It is also roughly orthogonal to the disk like structure 
regularly reported for the Milky Way satellite system and in close alignment with MSl's Giant Stellar Stream. 
A similar analysis of the asymmetry of the M3 1 satellite distribution finds that it is also significantly larger than 
one would expect from a random distribution. In particular, it is remarkable that 20 of the 27 satellites most 
likely lie on the Milky Way side of the galaxy, with the asymmetry being most pronounced within the satellite 
subset forming the aforementioned disk. This lopsidedness is all the more intriguing in light of the apparent 
orthogonality observed between the satellite disk structures of the Milky Way and M3 1 . 
Subject headings: gSilsLxics: distribution — galaxies: dwarf — galaxies: individual (M31) — galaxies: satellites 



1. INTRODUCTION 

The possibility that irregular distributions of satellite galax- 
ies may be a common feature of large galaxy halos was orig- 
inally bolstered by several studies of the anisotropic distri- 
bution of our own galaxy's satellites. Lynden-Bell (1976) 
found that the Magellanic Stream along with Sculptor and 
the Draco-Ursa Minor Stream and their associated dwarf 
spheroidal galaxies all appear to lie in the orbital plane of 
the Magellanic Clouds. In Lynden-Bell (1982), all the then 
known dwarf spheroidal companions of the Milky Way are 
identified as lying in one of two streams. Kroupa, Theis, & 
Boily (2005) examined the likelihood of producing the ob- 
served disk-like distribution of Milky Way satellites from a 
spherical or oblate dark matter halo. From comparisons with 
theoretical isotropic satellite distributions produced from such 
a halo, they find that the chance of producing the observed dis- 
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tribution from the dark-matter sub-halos of cold-dark-matter 
(CDM) cosmology is less than 0.5 %. They examine various 
combinations of the inner most satellites and find a best-fit 
plane that is almost perpendicular to the plane of the Milky 
Way with a root-mean- square height ranging from only about 
10 to 30 kpc. Zentner et al. (2005), whilst finding a similar 
plane to Kroupa, Theis, & Boily (2005) for the Milky Way 
satellites, disagree with their assumption that such a plane 
is unlikely to arise from a conventional CDM dark matter 
halo. They argue that that the most luminous satellites cannot 
be taken for granted as forming randomly from the isotropic 
sub-halo distribution but instead, lie preferentially at smaller 
distances from the halo centre and co-planar with the major 
axis of the host halo. Coupled with the finding that galaxies 
preferentially align themselves with their major-axis highly- 
inclined or even perpendicular to that of the surrounding mat- 
ter (e.g. Navarro, Abadi, & Steinmetz 2004; Hartwick 2000), 
this would provide a good explanation for the observed ori- 
entation of the best fit plane. The arguments of Zentner et al. 
(2005) are contended by Metz, Kroupa, & Jerjen (2007) how- 
ever, and it should be noted that in contrast to Zentner et al. 
(2005), Libeskind et al. (2005) found that the distribution of 
the most massive subhaloes is not as flattened as the distribu- 
tion of the Milky Way's satellites. 

More recently, Lovell et al. (201 1), using the six halo mod- 
els in the Aquarius Simulations (Springel et al. 2008), find 
that all six halos produce a significant population of sub-halos 
with quasi-planar orbits aligned with the main halo spin. This, 
they argue, is a natural explanation for the observed satellite 
distribution of the Milky Way. Pawlowski et al. (2012) ar- 
gue against this however. With the calculation of the angular 
momenta of 8 Milky Way Satellites (Metz, Kroupa, & Libe- 
skind 2008) revealing a strong alignment between 6 of the 
orbital poles, Pawlowski et al. (2012) examine the likelihood 
of randomly drawing 6 sub-halos from each of the 6 Aquar- 
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ius simulations (among other halo simulations), and finding a 
similar degree of alignment. More precisely, they draw 10^ 
sets of 8 satellites from each of the 6 simulations, and se- 
lect the 6 with the highest degree of alignment between their 
orbits, thus emulating the findings of Metz, Kroupa, & Libe- 
skind (2008). They then look at the degree of clumping of the 
orbital poles IS^sph as well as the angular distance of the aver- 
age of the orbital pole inclinations from the model equator d 
and find that the actual degree of planarity observed for the 
six satellites identified by Metz, Kroupa, & Libeskind (2008) 
(A^^ = 35.4° and duw = 9.4°) are equalled or exceeded in 
the random draws in less than 10% of cases when Agph is con- 
sidered and less than 15% of cases for d. Starkenburg et al. 
(2012) also find that the degree of planarity observed for the 
Milky Way satellites is uncommon in all six of the Aquarius 
halos (see Fig. 7 of that study). 

In addition to the revelation that the Milky Way's satellites 
appear to inhabit highly-inclined great planes, they also ap- 
pear to corroborate the finding of Holmberg (1969), namely 
that the companions of Spiral Galaxies preferentially congre- 
gate at high galactic latitudes (the Holmberg Eff'ect), as ob- 
served in his study of 174 galaxy groups. It is not clear why 
this should be the case, or even if it truly is the case, although 
if the apparent adherence of satellite systems to polar great 
planes is typical of galaxies in general, then the Holmberg 
Eff'ect seems to be an extension of this. Quinn & Goodman 
(1986) proposed that dynamical friction may be responsible 
for the observed polar great planes, with those orbits spend- 
ing the most time in close proximity to the galactic disk, ex- 
periencing the fastest decay, while those that take the most 
direct route through the disk environs, namely the polar or- 
bits, experiencing the slowest orbital decay. It should be noted 
however that dynamical friction, whilst producing more polar 
orbits, would not produce planes of satellites. Nor would it 
be eff'ective on the young globular clusters which are shown 
to be co-planar with the "Vast Plane of Satellites" identified 
around the Milky Way by Pawlowski, Pflamm-Altenburg, & 
Kroupa (2012B). Indeed, Angus et al. (2011) show that dy- 
namical friction did not play a role in the formation of the 
Milky Way satellite orbits. 

Besides the conjecture that satellite great planes trace the 
major- axis of the dark-matter halo in which the parent galaxy 
resides, there are other proposed mechanisms for their cre- 
ation. One hypothesis is that these planes trace the orbits 
of ancient galaxies that have been cannibalized by the host 
galaxy. Palma, Majewski, & Johnston (2002) have inves- 
tigated this hypothesis by looking for planes among groups 
of satellite galaxies and globular clusters in the Milky Way's 
outer halo and find various members to be co-planar with ei- 
ther the Magellanic or Sagittarius streams. The findings of 
Lynden-Bell & Lynden-Bell (1995) are also consistent with 
such a hypothesis. Indeed, it is this hypothesis which is most 
strongly supported by Pawlowski et al. (2012), wherein the 
Asph and d of satellites drawn from various tidal models equal 
or exceed AfJ^ and dMw in over 80% of draws in some cases. 
A similar hypothesis, which in some regards links the galaxy - 
cannibalization and dark-matter hypotheses, proposes that the 
observed planes result from the orientation of the large-scale 
filamentary structure of galaxy clusters (e.g. Knebe et al. 
2004), an orientation traced out by those minor galaxies which 
fall into the halo of a major galaxy. Metz et al. (2009) argue 
however that extra-galactic associations of dwarf galaxies are 
too extended to account for the high degree of planarity ob- 



served for the Milky Way satellites. This argument is sup- 
ported by the findings of Vera-Ciro et al. (201 1) based on the 
Aquarius Simulations. 

The great obstacle to a conclusive resolution of these issues 
is the lack of systems for which reliable spatial (and kine- 
matic) data exists. While some such data does exist for large 
galaxy clusters such as Virgo and Coma, accurate 3D distri- 
butions of galaxies within their halo have for a long time been 
known only for our own galaxy's halo, ascertainable due to 
our central position within it. It has only been in recent times 
that a second system has opened up to us - that of our coun- 
terpart in the Local Group, M31. Whilst various databases 
of photometry and other data have been available for M31 
and some of its brighter companions for over a decade, it is 
the Pan- Andromeda Archaeological Survey (PAndAS - Mc- 
Connachie et al. 2009) - a deep photometric, 2-colour survey 
providing a uniform coverage of the M3 1 halo out to approx- 
imately 150 kpc - that has provided a new level of detail for 
this system. It is from this survey that we obtained our dis- 
tances to M3 1 and 27 of its companions, following the method 
developed in Conn et al. (2011) (henceforth CLIll) and fur- 
ther adapted for this purpose in Conn et al. (2012) (henceforth 
CIL12). The distances themselves and their associated un- 
certainty distributions are presented in CIL12 and it is these 
distributions that are utilized for all analysis contained in this 
paper. 

With regard to previous studies of the anisotropy in the M3 1 
satellite distribution, two investigations warrant consideration 
at this point. McConnachie & Irwin (2006), making use of 
Wide Field Camera (WFC) photometry from the Isaac New- 
ton Telescope (INT) in what was essentially the forerunner 
to the PAndAS Survey, focus on "Ghostly Streams" of satel- 
lite galaxies following a similar approach as Lynden-Bell & 
Lynden-Bell (1995) used for the Milky Way. In addition, 
they characterize the large degree of asymmetry in the satel- 
lite distribution, a feature also noted in CIL12, and examine 
the radial distribution of the satellites, noting a (statistically 
insignificant) larger average distance from M31 than that ob- 
served between the Milky Way and its satellites. They find 
a large number of candidate satellite streams, with some fa- 
voring the dwarf spheroidal members. Koch & Grebel (2006) 
utilize distance measurements from a variety of sources and 
focus particularly on planes of satellites and, whilst they do 
not find a particularly significant best fit plane when their 
whole satellite sample is considered, it is rather interesting 
that they find a 99.7 % statistical significance to their best fit 
plane when the then-known dwarf spheroidal galaxies domi- 
nate their sample. Furthermore, this plane is near-polar - as 
has been observed for the Milky Way, although they find lit- 
tle support for the Holmberg Eff'ect. Koch & Grebel (2006) 
utilize a particularly robust method in their search for high- 
significance planar fits to subsets of galaxies by considering 
every possible combination of a given number of satellites 
from their sample. 

In the current study we employ a similar approach, but with 
the great advantage of having a considerably extended sample 
of galaxies in our sample, with all distances derived by the 
same method and from the same data as described in CLIll 
and CIL12. (Distances are sampled from the distributions 
without the density prior applied - e.g. CIL12, Table 1, Col- 
umn 2 - whilst the best fit distances are drawn from CIL12, 
Table 2, Column 4.) As a result, we are able to give full con- 
sideration to the eff'ects of selection bias on the observed satel- 
lite distribution. This then presents an excellent opportunity 
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to greatly improve our knowledge of the three-dimensional 
structure of the M3 1 satellite distribution, with important im- 
plications regarding the recent evolution of the system. 

A breakdown of the structure of the paper is as follows. In 
Section §2, we outline our method for plane fitting (§2.1) and 
locating significant planes of satellites as well as the orien- 
tation, magnitude and significance of the asymmetry of the 
distribution. A method for generating random realizations 
of satellites subject to the same selection biases as the real 
data is also discussed in this section (§2.2) as is the selec- 
tion bias itself (§2.3). §3 then presents the results of apply- 
ing these methods, first to the sample as a whole, and then to 
subsets of galaxies. Specifically, §3.1 presents a study of pla- 
narity within the satellite system when all satellites contribute 
to the determination of the best fit plane; §3.2 examines the 
asymmetry in a similar way; §3.3 examines the orientations 
of planes of smaller subsets of satellites within the distribu- 
tion; and §3.4 concludes this section with a determination of 
the significance of a 'Great Plane' of satellites emerging from 
the preceding sections. Sections 4 and 5 then follow with dis- 
cussion and conclusions. 

Note that this paper was written in conjunction with a 
shorter contribution (Ibata et al. 2013; hereafter ILC 13) which 
announced some of the key discoveries resulting from the 
analysis we present here. In particular, the process of iden- 
tifying the member satellites of the 'great plane' discussed in 
ILC 13 is described here in more detail. In this analysis how- 
ever, we concern ourselves with the spatial structure of the 
satellite system only and so the reader should refer to ILC 13 
for the interesting insight provided by the addition of the ve- 
locity information. 

2. METHOD 

2.1. Plane Fitting 

In order to find planes of satellites within the M3 1 satellite 
system, our first concern is to convert the satellite distances 
as presented in CIL12 into three-dimensional positions. To 
do this, we begin with an M31 -centered, cartesian coordinate 
system oriented such that the x and y axes lie in the M3 1 tan- 
gent plane with the z-axis pointed toward the Earth. Specifi- 
cally, the X-axis corresponds to T]tp = which is the projection 
of M3rs Declination onto the tangent plane. The y-axis then 
corresponds to ^tp = - the projection of M3rs Right As- 
cension onto the tangent plane. The z-axis then points along 
the Earth-to-M3 1 vector, with magnitude increasing with dis- 
tance from Earth. This orientation can be seen in Fig. 10(c) 
of CIL12. Thus: 

X = DsatC0s(6)tSin(^) 

y = Dsatsin(r]) (1) 

Z = DsatCOS(e) - Dm31 

where Dmsi and Dsat are the distances from Earth to M31 
and from Earth to the satellite respectively, 6 is the angular 
separation on the sky between M3 1 and the satellite, and rj and 
^ are the real-angle equivalents of the tangent plane projection 
angles rjtp and ^tp respectively. 

Next, we rotate this reference frame to the conventional 
M3 1 reference frame such that the positive z-axis points to- 
ward M3rs north galactic pole^ (i.e. Zpmsi = +90°) and the 
Im31 = 0° meridian passes through the Earth. So as to be con- 
sistent with the earlier work of McConnachie & Irwin (2006), 
we have adopted the same values for M3rs position angle 



(39.8°) and inclination (77.5° - de Vaucouleurs 1958). Each 
object is hence rotated by 39.8° about the z-axis to counter 
the eff'ect of its position angle, and then 77.5° about the x- 
axis to account for M3rs inclination. A final rotation of 90° 
about the z-axis is then necessary to bring /^si = 0° into 
alignment with the direction of Earth (which hence lies at 
Im31 = 0°,Z7m31 = -12.5°). The resulting spherical coordi- 
nates for each object in the sample are plotted onto an Aitoff'- 
Hammer projection in Fig. 1. This same figure also shows the 
uncertainties in position associated with each object, gener- 
ated via sampling of the respective distance posterior proba- 
bility distributions (PPDs) of each object and subsequent con- 
version of each drawn distance into a three-dimensional posi- 
tion. 

With the satellites' positions determined in cartesian coor- 
dinates, it is straight forward to determine the minimum dis- 
tance of each satellite from a given plane as follows: 

Dpiane = l^^ + by + cz + d\ (2) 

where D plane is the distance of a satellite at a point (x, y, z) 
from a plane whose normal vector is (a, c) and is of unit 
length. For simplicity, we invoke the reasonable requirement 
that all planes must pass through the center of M3 1 and so 
in our case, d = and the plane normal vector points out 
from the center of M3 1 . Hence, in order to find the best-fit 
or maximum significance plane to a set of satellites, we need 
simply minimize Dpiane for the satellites to be fitted. This can 
be done via a variety of means, some of which are compared 
in the following section, but perhaps the most robust and the 
predominant method employed in this study, is that of min- 
imizing the root-mean- square (RMS) of the distances to the 
fitted satellites. 

In order to measure the asymmetry of the satellite distribu- 
tion about a given plane, we need only count the number of 
satellites on one side of the plane. To do this, we can simply 
remove the absolute value signs from equation 2, so that the 
side of the plane on which a satellite lies can be determined by 
whether Dpiane is positive or negative. The plane of maximum 
asymmetry is then taken to be that which divides the sample 
such that the diff'erence in satellite counts for opposite sides 
of the plane is greatest. 

Whether we wish to determine the best fit plane through a 
sample of satellites or the plane of maximum asymmetry, we 
require a system by which a large number of planes can be 
tested on the sample so that the goodness of fit (or asymme- 
try) can be calculated for each. To do this, we define each 
tested plane by its normal vector or pole (a, b, c) so that Eq. 
2 can be applied directly. We then rotate this pole to diff'er- 
ent orientations around the sky in such a way as to 'scan' the 
whole sphere evenly and at a suitably high resolution. In prac- 
tice, we need to be able to apply this routine many thousands 
of times for a large number of samples and so a fast computa- 
tional time is of the essence. To this end, for a given sample, 
our algorithm determines the desired plane following a two 
step procedure. 

Firstly, a low resolution scan of the sphere is made to de- 
termine the approximate direction on the sky of the pole to 
the best-fit plane. Only half the sphere actually needs to be 
scanned since poles lying on the opposite hemisphere corre- 
spond to the identical planes flipped upside down. The low 
resolution scan tests 2233 difl'erent poles across the hemi- 

^ Defined so as to point north in Equatorial coordinates 
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Fig. 1. — An AitofF-Hammer Projection showing the positions of MSl's sateUites, along with their associated Icr uncertainties. The positions show where each 
object would appear in the sky if viewed from the centre of M31, and are given in MSl's inherent galactic latitude and longitude. The position of the Milky Way 
is also shown for reference. The position uncertainties trace single arcs across the M3 1 sky, rather than two-dimensional patches on account of the restriction of 
the uncertainty to lie solely along the line of sight to the object from Earth. For this same reason, all the lines point radially outward from the Milky Way. These 
uncertainties also take into account the uncertainty in MSl's distance. The existence of a prominent plane, broadly consisting of Andromedas I, XI, XII, XIII, 
XIV, XVI, XVII, XXV, XXVI, XXVII, XXX and NGC147 and NGC185, is hinted at by the close proximity of their respective arcs. 



sphere. A near-uniform coverage is achieved by decreasing 
the number of planes tested in proportion to the cosine of the 
latitude of the planes' pole. This prohibits what would oth- 
erwise be an increased coverage at the higher latitudes of the 
coordinate system. With the pole to the best-fit plane deter- 
mined in low-resolution, a high resolution search is then made 
around the identified coordinates at 10 times the resolution. 
In this way a pole can eff'ectively be found at any of approxi- 
mately 250,000 evenly spread locations on the hemisphere. 

2.2. Generating Random Satellite Samples 

Whilst we are now equipped to identify best-fit planes to 
our sample and subsamples thereof, it is necessary to have 
some means of determining the significance of these planes 
in an absolute sense. The most intuitive way to do this is to 
perform the same analysis on a randomly generated sample 
of equal size. In particular, when we are concerned with all 
possible combinations of a particular number of satellites that 
can be produced from the whole sample, we are often deal- 
ing with a very large number of subsamples and so it is in- 
evitable that some of these subsets of satellites will exhibit a 
very high degree of planarity. Identical analysis must there- 
fore be performed on random distributions, to see if there are 
similar numbers of subsets with equal degrees of planarity. 

For this reason, considerable care was taken to design an 
algorithm capable of providing a unique random realization 
of the desired number of satellites whenever it is called. The 
algorithm makes use of the distance PPD for each satellite, 
and also takes into account the irregular window function (i.e. 
useable portion) of the PAndAS survey. Each time a satellite 
is to be added to the random realization, one of the 27 actual 
satellites is chosen at random and a distance is drawn from 



its associated PPD. This distance (Dsat) is then converted into 
a three dimensional position (x, y, z) following equation set 1 
and this satellite-to-M3 1 separation vector is then spun around 
to a new, random location in the M3 1 sky. Note that for each 
random realization, a new value of Dm3i is similarly drawn 
from the M31 distance PPD. 

Once again, care must be taken in this step to ensure that the 
whole sphere is given equal weight, otherwise there is a higher 
likelihood for the artificial satellites to be positioned at high 
latitude. Again, this is remedied by weighting the likelihood 
by the cosine of the latitude. 

With the new, random location for the satellite chosen, it 
is then projected back onto the sky as it would appear from 
Earth and a check is made to ensure that it does indeed lie 
within the boundaries of the PAndAS survey area, and outside 
of the central ellipse (5° major axis, 2° minor axis - see Fig. 
10 (c) of CIL12) where the disk of M31 inhibits reliable mea- 
surements. If the satellite does not meet these requirements, it 
is rejected and the satellite drawing process is repeated until a 
suitable position is generated. By repeating this process until 
the desired number of satellites are produced, a new, random 
comparison sample is generated which gives full account to 
the constraints on the actual data. 

In order for the random satellite realizations to mimic the 
actual data most closely, it is necessary that each artificial 
satellite is represented not by just one point, but rather a string 
of points reflecting the uncertainty in the Earth-to-Object dis- 
tance. Hence once acceptable positions for each satellite are 
drawn as described above, the distance distributions for each 
object are sampled and projected to their equivalent positions 
along the line of sight about the initially placed point. For 
sections 3.1, 3.2 and 3.4 each artificial satellite's distance dis- 
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tribution is represented by 1000 points such that each plane- 
fitting measurement is made for 1000 possible positions of 
the object and then the average value of the measurements 
is taken. The only exception to this number is where the 
maximum-likelihood approach is used in §3.1. Due to the 
inclusion of a second fitting parameter in this case, only 100 
samples are taken for each satellite. For §3.3, as we are not 
concerned with comparisons of plane significance between 
the real sample and the random realizations, it is sufficient 
to use a single drawn position for each artificial satellite. 

2.3. A note on Satellite Detection Bias 

By employing a similar method to that described above, it is 
also possible to explore the effect of the PAndAS survey area 
boundaries on the satellite detection bias as viewed from the 
center of M3 1 . It is intuitive that more satellites are likely to 
be detected along the line of sight to Earth, since even satel- 
lites at a large distance from M3 1 will still appear within the 
survey boundaries if they lie along this line. We can visual- 
ize this effect by generating a large number of randomly dis- 
tributed satellites and plotting them on the M3 1 sky after first 
rejecting those satellites that would appear outside the survey 
area 'mask' if viewed from Earth. To do this, one million 
satellites were drawn from a spherically symmetric halo po- 
tential with density falling off as a function of the square of the 
distance from the halo center. Satellites were hence drawn at 
distances between and 700 kpc from M3 1 with equal prob- 
ability. The satellites were then projected onto the M31 tan- 
gent plane and those satellites lying outside the survey area or 
inside the M31 disk obstruction area were excised from the 
density map. The resulting anisotropy of the satellites on the 
M31 sky is presented in Fig. 2. 
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Fig. 2. — An Aitoff-Hammer projection illustrating the satellite detection 
bias resulting from the PAndAS survey boundaries and M31 disk obstruction. 
Note that this figure utilizes a Gaussian blurring of radius 5°, as do all of the 
subsequent pole-density plots. 



As can be seen from the figure, the probability of detection 
is indeed higher along a great circle oriented edge-on with re- 
spect to the direction of Earth, and perpendicular to the M3 1 
disk (Zpmsi = 0°)- This great circle has its pole/ anti-pole at 
/m3i = ±90°,Z7m31 = 0° and hence we would expect a pre- 
disposition toward finding planes of satellites with a pole in 
this vicinity. We would also expect, though to a lesser ex- 
tent, to find an excess of satellite planes oriented edge-on 
with respect to Earth at any inclination. Such planes would 
have poles lying anywhere on the great circle whose normal 
is directed toward Earth. The drop in the satellite density at 
/m31 = 0°,/?M3i = -12.5° and/M3i = ±180°,/?m31 = 12.5° 
is a consequence of the hinderance to detection caused by the 
M3 1 disk. Due to the increased volume of space covered by 



the survey at greater distances from Earth, unhindered satel- 
lite detection is possible over a larger range of angles on the 
far side of M31 in comparison to the Earth- ward side. 

3. RESULTS 
3.1. Best Fit Plane to the Entire Satellite Sample 

In order to find the best-fit plane to the satellite system as 
a whole, the procedure of §2.1 is applied to the whole sam- 
ple of 27 satellites presented in CIL12. The RMS thickness 
of the sample is used here, as in subsequent sections, as the 
statistic of planarity; we find it to be a robust measure and 
it has the convenient property of being computationally inex- 
pensive. Since we are dealing with only one sample in this 
case, two other measures are also used for comparison. The 
first calculates the sum of the absolute values of the distances 
of each of the satellites from the tested plane. The second is 
essentially a maximum likelihood approach and replaces the 
plane of zero-thickness with a 'Gaussian Plane' such that a 
satellite's position within the Gaussian determines the plane's 
goodness-of-fit to that satellite. This second approach re- 
quires that diff'erent Gaussian widths a be tested for each 
plane orientation in order to find the width that best matches 
the satellite distribution. Values between 5kpc and 150 kpc 
were tested at 5 kpc intervals for each tested plane orienta- 
tion. Hence an additional characteristic of the satellite distri- 
bution is obtained, but at the expense of a considerably longer 
computation time. 

For each of the three measures of goodness-of-fit described 
above, the first step is to find the best-fit plane to the satel- 
lite positions with their positions determined from their best- 
fit distances. When either the RMS or maximum likeli- 
hood approach is used, the same best-fit plane is found as 
0.153JC -h 0.932}; + 0.329z = with pole at (/m31,^m3i) = 
(-80.7°, 19.2°). This plane is plotted as a great circle on the 
M31 sky in Fig. 3 with the poles of the plane indicated. 
When the absolute distance sum is used instead, the pole is 
found farther from the plane of the galaxy, at (/m31, ^m3i) = 
(-74.9°, 24.3°). Nevertheless, the polar-plane described by 
Koch & Grebel (2006) is supported by either measurement, 
and is reminiscent of the satellite streams identified in the 
Milky Way satellite system. In light of the detection biases 
imposed by the PAndAS survey area as illustrated in Fig. 2, 
the result in this case must clearly be treated with suitable 
caution however. Like Koch & Grebel (2006), we find little 
evidence for the Holmberg Eff'ect, with only 3 best-fit satel- 
lite positions falling within 30° of the M31 galactic poles, 
and only 6 of the Icr error trails from Fig. 1 pass beyond 
Z?M31 = ±60°. 

To determine the uncertainty in the plane's goodness-of-fit, 
we need to repeat the procedure for a large number of real- 
izations of the satellite sample, with the best-fit satellite dis- 
tances replaced with a distance drawn at random from their re- 
spective satellite distance PPDs. A density map of the best-fit 
plane poles identified from 200,000 such realizations is pre- 
sented in Fig. 4. This figure was generated using the distribu- 
tion RMS as the goodness-of-fit statistic, and contains 71.1% 
of all poles within a 5° radius of the best-fit pole stated above. 
When the sum of absolute distances is used in place of the 
RMS, this fraction falls to 68.3%, or to 70.9% when the max- 
imum likelihood approach is used. It should be noted that 
the distribution of poles lies in close proximity to the pole of 
maximum detection bias at /m3i = -90°,/7m31 = 0°, again 
suggesting that the detection bias is having a strong influence 
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Fig. 3. — An AitofF-Hammer Projection showing the best-fit plane to the 
satellite system as a whole. The pole and anti-pole of the plane are denoted 
by and 'x' symbols respectively. Only the best-fit satellite positions were 
incorporated into the fit for this figure. The distribution of poles obtainable 
from other possible realizations of the satellite distribution is presented in Fig. 
4. Note that the plane is near-polar, similar to the preferred plane orientations 
identified for the Milky Way Satellite System. 



on the polar orientation of the best-fit plane. 




Fig. 4. — A pole-density map showing the effective uncertainty in the loca- 
tion of the best-fit plane to the whole satellite sample. The poles of the best-fit 
planes derived for 200,000 possible realizations of the data are plotted, along 
with their corresponding anti-poles. 



In order to determine whether the goodness-of-fit of the 
best-fit plane is really physically significant, similar analysis 
should be performed on a large number of random realiza- 
tions of satellites, to see how often distributions of satellites 
arise with a comparable degree of planarity. Figure 5 presents 
probability distributions of the plane significance for possi- 
ble realizations of the real satellite sample along with average 
values from random realizations of the satellites (as per §2.2), 
obtained using the three measures of goodness-of-fit stated 
above. 

It is immediately clear from Fig. 5 that regardless of the 
choice of the measure of goodness-of-fit, the range of values 
obtainable from possible realizations of the real satellite po- 
sitions are similar to the most likely values to be expected 
from completely random realizations of the satellites. Hence, 
whilst a prominent plane of satellites comprising roughly half 
of the sample is suggested in Fig. 1, it would seem that the 
sample as a whole is no more planar than would be expected 
from a strictly random distribution. Again, this is in keep- 
ing with the findings of Koch & Grebel (2006), and detracts 
from any physical significance that should be attributed to the 
plane's polar orientation. 

Further to this finding, the overall width of the 'plane' is 



again in keeping with that expected from a purely random 
satellite distribution. From fitting the Gaussian Plane to the 
best-fit satellite positions, a Icr width of 60 kpc is found to 
produce the best fit to the data. When the 200,000 PPD- 
sampled realizations were tested, a 1 cr of 60 kpc was found 
preferential in 66.3% of cases, with a 1 cr of 55 kpc being pre- 
ferred in 32.7% of cases. Values of 50 kpc make up the re- 
maining 1 % almost entirely. The average value for the actual 
satellite distribution was thus determined as 58.3 kpc. This 
value is similar to the most likely width identified from the 
10,000 random realizations, as can be seen in Fig. 6. 

3.2. The Plane of Maximum Asymmetry 

To determine the plane of maximum asymmetry and its sig- 
nificance, we employ an identical approach as in the preced- 
ing section, but with the goodness-of-fit statistic replaced with 
a count of the number of satellites on each side of the plane 
as per §2.1. As was suggested by the three-dimensional satel- 
lite distribution generated in GIL 12, the asymmetry about the 
M31 tangent plane is close to a maximum, with 19 satellites 
on the near- side of the plane but only 8 on the other when the 
best-fit satellite positions are assumed. The highest asymme- 
try plane possible from this same distribution has 21 satellites 
on one side and 6 on the other, with the equation of the plane 
identified by the algorithm as -0.797x- 0.315};-h 0.515z = 0. 
The anti-pole of this plane lies 27.2° away from the Milky 
Way at (/m3i, ^m3i) = (-21.6°, -31.0°). This plane is plotted 
as a great circle on the M31 sky in Fig. 7. 

When 200, 000 realizations of the satellite sample are gen- 
erated using the satellite's respective distance probability dis- 
tributions, the most likely asymmetry of the sample is actually 
found to be greater than this, with 23 satellites on one side 
and only 4 on the other. Such a scenario is more than twice 
as likely as the 21 : 6 scenario. In one realization, a plane 
was identified which could divide the sample such that all 27 
satellites lay in a single hemisphere, while an asymmetry of 
26 : 1 was found possible for 815 (0.4%) of the realizations. 
The distribution of maximum- asymmetry poles on the sky, as 
determined from realizations of possible satellite positions, is 
illustrated in Fig. 8, whilst Fig. 9 {a) plots the probability 
distribution for the greatest number of satellites that can be 
found in one hemisphere for a given realization of the ob- 
served satellite sample. The average value of this distribution 
is 22.7 (shown as a dashed line in Fig. 9 (/?)), a value which 
is equalled or exceeded for 422 out of the 10, 000 random re- 
alizations represented in Fig. 9 ib). A maximum asymmetry 
ratio of 21 : 6, as was observed for the best-fit satellite dis- 
tribution plotted in Fig. 7, is more common however, falling 
inside the Icr credibility interval. 

What is particularly striking about the satellite distribution 
however, is the orientation of the asymmetry, with the ma- 
jority of satellites lying on the near-side of the M31 tangent 
plane. From Fig. 9 (c), it is clear that the eff'ect of the dis- 
tance uncertainties lying along the line of sight is to create 
quite a broad distribution in the level of asymmetry about the 
tangent plane, though the average is markedly high at 20.3. 
To investigate the likelihood of this scenario arising from a 
random satellite distribution, we measure the average num- 
ber of satellites on either side of the M3 1 tangent plane for 
each of 10, 000 random realizations as per §2.2. The results 
are illustrated in Fig. 9 d). The observed profile is more- 
or-less as expected, with a maximum probability close to the 
minimum possible asymmetry at 14 and then a rapid fall off 
toward higher asymmetries. It is therefore clear that the dis- 
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Fig. 5. — Probability distributions for the planarity of the entire sateUite sample, as determined from three different measures of the plane goodness-of-fit. The 
left-hand column of figures gives the distribution of the goodness-of-fit statistic as obtained via plane fitting to 200,000 separate samplings of the real satellite 
sample. The right-hand column of figures summarizes the same procedure performed for 1 , 000 separate samplings of each of 10, 000 random realizations of the 
satellites (as per §2.2). It is important to note that each histogram in this column has been generated by plotting the average values from the 10, 000 individual 
histograms corresponding to each of the random realizations and hence they should only be compared with the average of the histograms in the left-hand column. 
The goodness-of-fit statistic for a) and b) is the distribution RMS; for c) and d) is the absolute distance sum and; for e) and f) is the sum of satellite likelihoods. 
The average of the histograms in (a), (c) and (e) are shown in (b), (d) and (f) respectively as dashed lines. Red, green and blue lines denote the extent of Icr 
(68.2%), 90% and 99% credibility intervals respectively. 



tance uncertainties lying along the line of sight have no sig- 
nificant bearing on the orientation of the asymmetry. Yet the 
observed degree of asymmetry about the M3 1 tangent plane is 
equalled or exceeded in only 46 of the 10, 000 random satellite 
realizations and hence is very significant. The possibility that 
this asymmetry may be a consequence of data incomplete- 
ness is currently being examined more closely (see Martin et 
al. 2013), although it seems very unlikely. The high degree 
of asymmetry is still observed even when only the brightest 
satellites are considered. Furthermore, the data incomplete- 
ness appears to be dominated by the boundaries of the PAn- 
dAS survey area and obstructed regions which are already 
taken into account by our analysis. Indeed, one would ex- 
pect more satellites to be observed on the far side of the M3 1 
tangent plane on account of the increased volume of space 
covered by the survey at greater distances, an eff'ect clearly 



visible in Fig. 2. 

3.3. Subsets of Satellites 

It is perhaps not surprising that the satellite system of M3 1 , 
when treated as a whole, is no more planar than one would 
expect from a random sample of comparable size. Indeed, 
a similar result was noted for the M31 system by Koch & 
Grebel (2006). The existence of outliers in our satellite sam- 
ple was already clear from Fig. 1 and furthermore, if multiple 
planes of diff'ering orientation are present as has been sug- 
gested for both the Milky Way's satellite system (e.g. Lynden- 
Bell 1982; Pawlowski, Pflamm-Altenburg, & Kroupa 2012B) 
and the M31 system (McConnachie & Irwin 2006), then the 
goodness of fit of the best-fit plane to the entire distribution is 
of little consequence. For this reason, we now concentrate our 
analysis on subsets or combinations of satellites. Specifically, 
we perform a pole-count analysis by determining the pole of 
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the best-fit plane to every possible satellite combination of a 
particular size that can be drawn from the entire sample. 

A pole-count analysis is an excellent way of mapping the 
degree of prominence of various planes that exist within the 
distribution as a whole, whatever their orientation may be. 
The choice of combination size is not trivial however. The 
number of combinations s of 3. particular number of satellites 
k that can be drawn from the entire sample of n satellites can 
be determined as follows: 



(3) 



Plane Sigma (kpc) 



Fig. 6. — The probability distribution for the average Icr width as deter- 
mined from 10,000 random distributions of 27 satellites. This figure is gen- 
erated from the same run as Fig. 5 f) and is the result of marginalizing over 
the plane-orientation model parameters. 
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Fig. 7. — An Aitoff'-Hammer projection showing the plane of maximum 
asymmetry identified from the full sample of best-fit satellite positions. It 
divides the distribution such that 21 satellites lie in one hemisphere, but only 
6 in the other. The anti-pole of the maximum asymmetry plane lies just 28.1° 
from the Milky Way as viewed from the center of M3 1 . 




Fig. 8. — A pole-density map showing the eff'ective uncertainty in the lo- 
cation of the maximum asymmetry plane to the whole satellite sample. The 
poles of the maximum asymmetry planes derived for 200,000 possible real- 
izations of the data are plotted, along with their corresponding anti-poles. The 
elongated distributions that run through the pole and anti-pole determined 
from the best-fit distribution (see Fig. 7) arise due to the orientation of the 
uncertainty trails of the individual satellite positions, as presented in Fig. 1. 
Note that the probability of the anti-pole of the asymmetry lying within a 
couple of degrees of the direction of the Milky Way is close to a maximum. 



kl(n-k)l 



For reasons that shall be discussed shortly, we will eff'ectively 
be working with a sample of 25 satellite positions. It is clear 
from this equation however that with 25 satellites forming the 
entire sample, the total number of combinations that can be 
drawn may be very large, depending on the number of satel- 
lites forming the combinations. For instance, if ^ = 25 and 
k = 13, there are over 5.2 million possible combinations that 
can be drawn. Additionally, if we are to properly account for 
the uncertainties in the satellite positions, it will be necessary 
to sample from the distance distributions of each satellite a 
large number of times for every combination. Given that we 
must test every possible plane orientation (as per §2.1) for 
every rendition of every combination, the computation times 
can become impracticable. It is therefore necessary to limit 
our combination sizes as much as possible. We note however, 
that the final pole-plot distribution showing the poles of the 
best-fit planes to each combination, is not so dependent on 
the combination size as might at first be thought. 

With all the planes tested as per §2.1 having to pass through 
the center of M3 1 , the minimum number of satellites that can 
not be fitted exactly is 3. This is therefore the smallest com- 
bination size we consider. There are 2, 300 combinations of 3 
satellites that can be drawn from the full sample of 25 satel- 
lites. If we increase the combination size considerably to 7 
satellites, there are 480, 700 satellite combinations that can be 
drawn. Due to an excessive number of combinations beyond 
this point, this is the largest combination size we consider. 
But it is critical to note that even if we produce our pole-plot 
map from combinations of only 3 satellites we do not exclu- 
sively find planes consisting of 3 satellites. If a plane of 7 
satellites exists for instance, then by Eq. 3, such a plane will 
produce 35 poles at the same location on the pole plot, where 
a plane consisting of only 3 satellites would contribute only 
one pole. Conversely if we take combinations of 7 satellites, 
despite the larger number of possible combinations in total, 
we become less sensitive to planes made up of less than 7 
satellites. So in a sense, the combination size we choose de- 
pends on the satellite planes we wish to be most sensitive to. 
In practice, we have found that the smaller combination sizes 
of 3 and 4 satellites are particularly useful for identifying the 
lowest RMS planes congregating around the band of satellites 
visible in Fig. 1. The larger combination sizes of 5, 6 and 
7 satellites gradually shift toward finding planes closer to the 
best-fit plane to the entire satellite sample illustrated in Fig. 3. 

Noting these points, we proceed as follows. First, the num- 
ber of satellites per combination k is chosen (3 < k < 1) and 
then for each combination, distances are drawn for each of the 
satellites from their respective posterior distance distributions 
as provided in CIL12. To give a satisfactory representation 
of the form of the distributions, each combination is sampled 
100 times. As such, each satellite combination contributes not 
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Fig. 9. — Asymmetry probability distributions. The top two histograms plot probability distributions for the greatest number of satellites that can be found 
in one hemisphere, as generated from (a) 200,000 samplings of satellite positions possible from the data and; (b) the average of 1000 samplings from each of 
10,000 random realizations of the satellites generated as per §2.2. Figures (c) and (d) give the equivalent distributions when the maximum asymmetry plane is 
replaced with the fixed M31 tangent plane. As for Fig. 5, the histograms in the right-hand column should only be compared with the average of the corresponding 
histogram in the left column. The average value of the histograms of (a) and (c) are shown in (b) and (d) respectively as a dashed line. 



1 pole to the pole density map for the chosen combination size 
but 100, with the spread of poles relating the possible orienta- 
tions of the best-fit plane to the combination, given the error 
in the individual satellite positions. The contribution of each 
pole to the density map is also weighted by dividing it by the 
RMS of the best-fit plane it represents. Thus each pole does 
not contribute 1 count, but rather some fraction, depending 
on how good a fit the plane it represents is to the satellites in 
the combination. This fraction is also further divided by 100, 
since it represents only 1 % of the samples for the combina- 
tion, as just discussed. 

As stated above, it should also be noted that we eff'ectively 
limit the total number of satellites in our sample to 25 for all 
analysis in this subsection. This is to account for the bound 
group of satellites consisting of NGC147, NGC185 and And 
XXX (henceforth the NGC147 group - see Irwin et al. 2013). 
Since we suspect that these satellites orbit M3 1 as a group and 
since they all lie along the apparent plane identified in Fig. 1, 
it is preferable to treat the group as a single object when we 
are not concerned with measurements of the significance of 
particular planes. To do this, we take the luminosity weighted 
centre as an approximation for the center of mass, and treat 
this determined position as though it were the location of a 
single satellite. To calculate the luminosity weighted center, 
we can ignore the contribution from And XXX since it is neg- 
ligible compared with the contributions of the two dwarf ellip- 
ticals. From the Third Reference Catalogue of Bright Galax- 
ies (de Vaucouleurs et al. 1991), NGC185 is 0.2 magnitudes 
brighter than NGC147 in the V-band. Each time the NGC147 
group is chosen as one of the 'satellites' for a combination, 
distances to each of NGC147 and NGC185 are sampled from 
their respective distributions and the luminosity weighted cen- 



ter of the group is determined. As for any other combination, 
this position, along with all other satellites in the combination, 
is sampled 100 times. 

The results of applying the above procedure to all combi- 
nations of 3, 4, 5, 6 and 7 satellites that can be drawn from 
the total sample is presented in Fig. 10. The left-hand column 
shows the fit to the most planar combination determined from 
the best-fit positions whilst the right-hand column shows the 
corresponding pole density plots for all combinations of that 
particular number of satellites, based on 100 samples of each 
combination as per the discussion above. It is noteworthy that 
the best-fit planes to the most planar combinations are almost 
identical in every case, except for that of the 3 satellite com- 
binations, where the RMS values are so small for so many 
combinations as to make this result not particularly impor- 
tant. It should also be noted that these best-fit planes trace 
out the same approximate great circle as the prominent plane 
indicated in Fig. 1, a result that shall be investigated a little 
later in §3.4. It is particularly interesting that the pole shared 
by each of these planes, located at /m3i ~ -80°,/?m31 ~ 40° 
corresponds to a pole count maximum in each of the pole 
plots. This indicates that many of the satellite combinations 
are aligned along this plane, hence further suggesting that the 
plane applies to more satellites than the combination sizes 
tested here. The other, lower latitude principle maximum in 
the pole plots is that corresponding approximately to the best 
fit to all the satellites and hence it grows more prominent in 
the plots made from larger combination sizes as discussed ear- 
lier. 

Besides the pole count maxima that are strongly indicative 
of a highly planar subset of satellites, the other principle fea- 
ture of the pole plots in Fig. 10 is the great circle along which 
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the pole count density is highest. This great circle is very 
prominent but great caution must be exercised in attributing 
any significance to it. It is centered on the Milky Way in- 
dicating that the constituent poles result from a majority of 
satellites lying along the Earth to M31 line of site. But this 
reflects the anisotropy predicted from Fig. 2, the result of the 
bias incurred by the finite area of the PAndAS survey. Hence 
it would seem that the progenitor of this prominent great cir- 
cle is not physical but rather the result of selection eff'ects. To 
further investigate the significance of the patterns observed in 
the pole plots, 1000 random realizations of 25 satellites were 
generated as per §2.2, and a similar pole count analysis per- 
formed on each of them. Specifically, the pole density distri- 
bution resulting from the best fit planes to all combinations 
of 5 satellites was generated for each of them. The resulting 
pole plots for 8 of the 1000 random realizations (chosen at 
random) are presented in Fig. 11 along with an enlarged ver- 
sion of the equivalent plot from Fig. 10 generated from the 
real distribution. A bias toward a similar high-density great 
circle is indeed observed in these plots, but the plot generated 
from the actual data features a conspicuously narrower great 
circle, and a much more constrained distribution in general. 
This appears to be primarily the result of the large fraction of 
satellites that lie along the prominent plane that is repeatedly 
identified and plotted in the left-hand column of Fig. 10. It 
should also be noted that this plane, whilst being oriented per- 
fectly edge-on with respect to the Earth, contains a significant 
fraction of satellites lying well outside the region of the M31 
sky where the detection bias is large, and hence it is unlikely 
that its prominence is due to our observational constraints. 

Figure 12 provides for a comparison between the concen- 
tration of poles around the principle maximum in the pole dis- 
tributions of the actual satellite distribution and the average of 
the 1000 random satellite distributions. From line (a) in Fig. 
12 we see that 21.5% of all combinations of the actual satel- 
lite positions are fitted by a best-fit plane with pole within 15° 
of the principal maximum (located at /m3i = -78.7°, Zpmsi = 
38.4°). This is in stark contrast to the 12.0% that lie within 
15° of the principal maximum for the average random real- 
ization of satellite positions (Fig. 12 line (b)). Furthermore, 
we find that only 1 17 of the 1000 random realizations exhib- 
ited the degree of concentration of poles within 15° of the 
principal maximum that was observed for the actual satellite 
distribution. Hence it would seem that a large percentage of 
satellite combinations are fitted by best-fit planes that all have 
strikingly similar orientations when compared with what one 
could expect from a random distribution of satellites. Again, 
this points toward a significant plane of satellites that includes 
a large fraction of the whole satellite sample. 

In order to obtain a better understanding of the satellites that 
this plane consists of, it is of particular interest to explore the 
number of times each satellite is included in a combination 
that is best fit by a plane with pole in close proximity to the 
principal maximum in the pole distribution for the entire sam- 
ple. Once again, we use the pole distribution for all combina- 
tions of 5 satellites, and we count the number of times each 
satellite contributes to a pole within 3° of the principal maxi- 
mum at /m3i = -78.7°, /?M3i = 38.4°. The counts are divided 
by 100 to account for the 100 samples that are taken of each 
combination. The result can be seen in Fig. 13. From this fig- 
ure, it can be seen that the main contributors to the principal 
maximum in pole counts are those same satellites identified 
as forming a prominent plane in Fig 1 , namely Andromedas I, 
XI, XII, XIII, XIV, XVI, XVII, XXV, XXVI, XXVII and the 



NGC147 group, along with Andromeda III and Andromeda 
IX. Hence the conclusion of our analysis thus far must be 
that there is indeed a significant plane in the satellite distri- 
bution of M3 1 and that it broadly consists of the aforesaid 
satellites. We therefore investigate the numerical significance 
of the best-fit plane to these satellites in §3.4. As yet there 
is still more to be gleaned from a study of the pole density 
distribution however. 

From Fig. 13 we have been able to determine the principle 
contributing satellites to the principal maximum in the pole 
density distribution, but what of the remaining satellites? Do 
the positions of these satellites follow any particular trend? 
The best way to determine this is to construct pole density 
plots of the two halves of the complete sample, namely the 
major contributors to the principal maximum and the minor 
contributors. The resulting pole plots are presented in Fig. 
14. 

The left-hand plot of Fig. 14 shows the pole density dis- 
tribution generated from the major contributing satellites to 
the principal maximum at /m3i = -78.7°,Z?m31 = 38.4°. 
This half-sample includes Andromedas I, III, IX, XI, XII, 
XIII, XIV, XVI, XVII, XXV, XXVI, XXVII and the NGC147 
group. As expected, this plot reflects the existence of the 
aforementioned plane with all combination poles lying in 
the vicinity of the principal maximum. The right-hand plot, 
with poles generated from the remaining 12 satellites, namely 
Andromedas II, V, X, XV, XVIII, XIX, XX, XXI, XXII, 
XXIII, XXIV and M33, paints a very difl'erent picture how- 
ever. There is a much greater spread in the distribution of 
poles, with the great circle induced by the survey area bias 
once again conspicuous. Also prominent in this plot are 2 
density maxima with their corresponding mirror images in the 
opposite hemisphere. The maximum lying midway between 
Andromedas XIX and XX lies very close to the pole of maxi- 
mum detection bias at /m3i = -90°, Bmsi = 0° and so it is not 
unexpected, now that the prominent plane of satellites is ef- 
fectively removed from the distribution. The elongated max- 
imum passing through /m3i ~ 45°, bM3i ~ 45° is more inter- 
esting however, and suggests the possibility of a second plane, 
roughly orthogonal to the major plane represented in the left- 
hand plot, though much less conspicuous. The planes repre- 
sented by this maximum pass close to the error trails on the 
M31 sky of Andromedas II, III, XIX, XX, XXIII and XXIV. 
This maximum is faintly discernible in the pole distribution 
for combinations of 6 satellites presented in Fig. 10 but is no 
more pronounced than anywhere else along the high-density 
great circle in any of the other pole plots. On account of this, it 
would appear that this plane is likely no more significant than 
one would expect to find from a random satellite distribution 
subject to the same detection biases, such as those illustrated 
in Fig. 11. 

3.4. A Great Plane of Satellites 

Throughout the investigation undertaken thus far, all evi- 
dence has repeatedly pointed toward a conspicuously planar 
sub- set of satellites consisting of roughly half the total sam- 
ple of satellites. Andromedas I, XI, XII, XIII, XIV, XVI, 
XVII, XXV, XXVI, XXVII and XXX as well as the dwarf 
ellipticals NGC147 and NGC185 all appeared to He along a 
plane in Fig. 1. The reality of this co-planarity was veri- 
fied in §3.3 and in particular Fig. 13, which also suggested 
that Andromeda III and Andromeda IX should be considered 
as plane members. Hence it is of great interest to ascertain 
whether this 'great plane' is in fact significant. To do this, it 
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Fig. 10. — Best fit planes and pole density maps for combinations of 3 through 7 satellites. The left-hand column shows the best-fit plane through the combination 
of satellites that can be fit with the lowest RMS. Satellites included in the best-fit combination are colored red. The centre of the NGC147 group is marked with 
a circle, and lies on the best-fit plane in every case. The three members of this group are colored orange. Only the best-fit satellite positions are considered for 
these plots. The right-hand column shows the corresponding pole density plot for the poles of all satellite combinations. These plots have been weighted by the 
RMS of each pole and fully account for the uncertainty in the satellite positions. 




Fig. 11. — Pole density maps for 8 random realizations of 25 satellites. The maps plot the poles for the best-fit planes to all combinations of 5 satellites. The 
contribution of each pole is weighted by the RMS of the plane it represents. The map resulting from all combinations of 5 satellites drawn from the real data 
is shown again at the top for comparison. Great circle overdensities are evident to varying degrees in the plots and are a result of the satellite detection bias 



The Distribution of the M3 1 Satellite System 



13 
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Fig. 12. — Radial density profiles showing the percentage of all poles ly- 
ing within n degrees of the densest point in the pole count distributions (the 
principal maximum) for a) the actual satellite distribution and b) the aver- 
age of 1000 random satellite distributions. The profile for the actual satellite 
distribution is generated from the same pole distribution as illustrated for 5 
satellites in Fig. 10 and at the top of Fig. 11. Note that the relative linearity 
of (b) compared with (a) is simply a result of the averaging of a large number 
of individual profiles undertaken to produce the former. 
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Fig. 13. — Histogram showing the relative contribution of each satel- 
lite to the pole density within 3° of the principal maximum at /m3i = 
-78.7°, Z?M3i = 38.4°. The histogram is generated from the same pole distri- 
bution as illustrated for 5 satellites in Fig. 10 and at the top of Fig. 11. 



is necessary to determine how likely such a plane is to arise 
from a random satellite distribution subject to the same selec- 
tion biases. The plane itself and the satellites of which it is 
constituted are illustrated in Fig. 15. The plane shown is that 
calculated from the best-fit satellite positions and has equa- 
tion of the form: 0.158x + 0.769}; + 0.620z = with pole at 
(/m3i,^m3i) = (-78.4°, 38.3°). Note that for this section, we 
re-instate NGC147, NGC185 and Andromeda XXX as sepa- 
rate objects since we are again concerned with measurements 
of the significance of the planarity of the distribution. Our 
'great plane' thus consists of 15 satellites out of the entire 
sample of 27. 

Using the method of §2.2, we again generate 10, 000 inde- 
pendent random realizations of 27 satellites and seek the most 
planar combination of 15 satellites from each. For each ran- 
dom realization, we sample 1000 possible positions for each 



satellite as in previous sections and take the average value for 
the RMS of the best fit plane through the most planar com- 
bination. Since there are more than 17 million ways that 15 
satellites can be drawn from 27, and since we are not con- 
cerned with the orientation of each fitted plane as we have 
been in all previous sections, we depart from the plane fit- 
ting method of §2.1 for this section and instead proceed as 
follows. For each sample of satellite positions from each re- 
alization, 10, 000 randomized planes are generated and the 15 
closest satellites of the 27 to the plane are stored in each case 
and the RMS recorded. The lowest RMS achieved is hence 
taken to be that for the most planar combination of 15 satel- 
lites in the sample. These minimum RMS values from each of 
the 1000 samples of the particular random realization are then 
averaged to provide the best representation for the realization, 
given the distance uncertainties. Fig. 16 provides probabil- 
ity distributions in the RMS for the observed 'great plane' {a) 
together with those for the average RMS for the most planar 
combination from each random realization (b). The average 
RMS for the observed plane is plotted in (b) for comparison. 

As can be seen from Fig. 16, the RMS for the observed 
plane is very low compared to what one could reasonably ex- 
pect from a chance alignment. Indeed, the average RMS of 
12.58 kpc for the observed plane is found to be equalled or 
exceeded in only 36 out of the 10,000 random realizations. 
The chances of obtaining such a planar group of 15 satellites 
from a sample of 27 at random is thus estimated as 0.36%. 
Hence we can conclude from this test that the observed plane 
is very unlikely to be a chance alignment, but rather the re- 
sult of some underlying physical mechanism. Note that an in- 
dependent but equivalent investigation is presented in ILC13 
where such an alignment is found to occur in only 0.15% of 
instances. This is due to the larger central obstruction adopted 
in that analysis (19.6 vs. 7.9 sq. deg.) which rejects more 
satellites in close proximity to the plane pivot point (M31) 
where small plane distances are most likely. 

4. DISCUSSION 

Throughout the analysis conducted in §3, the presence of 
a prominent plane of satellites has been a consistent feature. 
This is not the first time that a significant plane of satellites 
has been identified from among the denizens of the M3 1 halo 
however. As early as 1995, Fusi Pecci et al. (1995) had noted 
that the majority of the then known satellites lay conspicu- 
ously in a plane oriented edge-on with respect to the Milky 
Way. Koch & Grebel (2006) identified a highly significant 
plane lying within 5° to 7° of being polar. Furthermore, they 
identify a subset of 9 satellites from this plane lying within 
a thin disk with an RMS of 16 kpc. Metz, Kroupa, & Jer- 
jen (2007) and later Metz, Kroupa, & Jerjen (2009) simi- 
larly identify a disk of satellites, this time not so markedly 
polar, with pole (in our coordinate system) at (/m31,^m3i) = 
(-70.2°, 32.9°). They find this disk to have an RMS height of 
39.2 kpc. This disk is clearly the same structure that we iden- 
tify here, being tilted by only 8.6° with respect to our 'great 
plane.' Our plane is found to have a much smaller RMS of 
just 12.34^Q 43 kpc however, despite including a comparable 
number of satellites. It is particularly noteworthy however, 
that their satellite sample is significantly diff'erent to that used 
here, with their disk including M32, NGC205, ICIO, LGS3 
and IC1613 - all of which lie outside the portion of the R\n- 
dAS survey region used in this study (see Fig. 10 (c) of 
GIL 12). Indeed, it is clear from Fig. 4 of McConnachie & 
Irwin (2006) that the galaxies M32, ICIO, LGS3 and IC1613 
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Fig. 14. — Pole density distributions generated from all combinations of 5 satellites possible from: Left) the satellites contributing significantly to the principal 
maximum at /m3i = -78.7°, bM3\ = 38.4° as per Fig. 13 and Right) the remaining 12 satellites. 




Fig. 15. — A Great Plane of Satellites consisting of Andromedas I, III, IX, 
XI, XII, XIII, XIV, XVI, XVII, XXV, XXVI, XXVII, XXX, NGC147 and 
NGC185. The plane shown is that derived from the best-fit satellite positions. 
The pole is located at (/m3i, ^m3i) = (-78.4°, 38.3°). 



all lie along the same great circle as our 'great plane' in Fig. 
15, as do their entire error trails. Their conformity along with 
Andromeda I to a thin disk is noted in the said paper as one of 
8 possible 'streams of satellites,' thus providing another early 
detection of the plane identified by this study. Both Majewski 
et al. (2007) and Irwin et al. (2008) also draw attention to the 
linear distribution of many of the plane-member satellites on 
the sky, a consequence of the edge-on orientation of the plane 
as indicated by the present study. The plane of Metz, Kroupa, 
& Jerjen (2009) does however include a significant number of 
satellites that, whilst included in our sample, we exclude due 
to their looser association with our plane. This then accounts 
for the much smaller RMS height observed in our study. 

Unlike previous studies of the M31 satellite system, we 
have a significant advantage in this study on account of the 
greatly improved sample of satellites available to us. Our 
sample is not only more numerous, but the positions are all 
determined via the same method applied to the same data 
as per CLIll and CIL12. We are thus aff'orded unprece- 
dented knowledge of the satellite detection biases, as well 
as the uncertainties in the object positions and have factored 
this knowledge into the analysis. An understanding of this 
bias is of particular importance when it comes to ascertaining 
the significance of any substructure identified, since a phys- 
ically homogeneous satellite distribution will inevitably ap- 
pear anisotropic after 'folding in' the selection function and it 
is important that we do not attribute physical significance to 
this anisotropy. 

Even after taking these eff'ects into account however, there 




Fig. 16. — Determining the significance of the observed 'great plane' of 
satellites (see Fig. 15). Figure (a) gives the distribution of possible values 
of the RMS obtainable from 200,000 realizations of possible positions of the 
15 plane members, given their respective distance probability distributions. 
Figure (b) plots the average RMS of the best fit plane through the most planar 
combination of 15 satellites for each of 10,000 random realizations of 27 
satellites. These satellites are subject to the same selection biases as the real 
data. As for Fig. 5, histogram (b) should only be compared with the average 
of histogram (a), which is plotted in (b) as a dashed line. It is thus clear 
that the planarity observed for our 'great plane' of satellites is very unlikely 
to arise by chance. The Icr (68.2%), 90% and 99% credibility intervals are 
shown as red, green and blue lines respectively. 



The Distribution of the M3 1 Satellite System 
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can be little doubt that the plane described in §3.4 is a real 
physical object. The component satellites extend well into the 
regions of low detection bias in Fig. 2 and the analysis of the 
last section makes it clear that such a thin disk of satellites has 
very little chance of arising within a random satellite distribu- 
tion of the same size, even when subject to the same obser- 
vation biases. Furthermore, it should be noted that the study 
of the plane's significance in §3.4 is likely to be conserva- 
tive, given that if the satellites M32, ICIO, LGS3, IC1613 and 
NGC205 were to be included in the analysis, the significance 
of our observed plane would likely grow still further. What 
is also particularly interesting is that subsequent research has 
shown 13 of the 15 objects to be co-rotating. This result is 
discussed in more detail in ILC13. 

What then could be the progenitor of this 'great plane'? 
The polar orientation one might expect to arise had the satel- 
lites formed within the dark matter halo or had the dynamical 
friction proposed by Quinn & Goodman (1986) had sufficient 
time to take eff'ect is not observed. Similarly, the findings 
of Metz et al. (2009) seemingly preclude the possibility that 
the structure might be the result of the accretion of an exter- 
nal galactic association. Furthermore, there is apparently no 
marked distinction in the metallicities of the disk members 
compared with the non-disk members as one might expect 
from this scenario, though it is possible that some of the non- 
disk satellites may share the same origin as the disk members. 
There remains however the possibility that the satellites trace 
out the tidal debris of a galaxy merger. This is a particularly 
interesting possibility, especially since the plane, when pro- 
jected onto the M31 tangent plane, is in close alignment with 
the Giant Stellar Stream. Indeed, Hammer et al. (2010) show 
that the Giant Stellar Stream could feasibly be the product of 
a major merger event that began around 9 Gyr ago, sustained 
by the returning stars from a tidal tail oriented similarly to our 
'great plane.' 

If a link is to be established, the observed asymmetry of the 
system must also be compatible with the tidal scenario how- 
ever. It is of particular interest that, of the 13 co-rotating satel- 
lites in the plane, all but one lie on the near side of the M3 1 
tangent plane. Indeed, if we removed all of the plane member- 
satellites from the system, the remaining satellite distribution 
would no longer be significantly asymmetric. With almost all 
of the satellites currently on the near side of M3 1 , it would 
seem on first consideration that the progenitor event could not 
have occurred substantially more than a typical orbital time 
ago or else the satellites would have had sufficient time to dis- 
perse. This suggests the event responsible must have occurred 
within the last 5 Gyr. Other studies however have supported 
the proposition that a group of tidal dwarf galaxies could sur- 
vive for extended periods whilst retaining the asymmetry in- 
herent from the time of their formation. For instance. Due 
et al. (2011) find three moderately old tidal dwarf candidates 
with this asymmetry preserved. It must also be remembered 
that the precise orbits of the satellites are undetermined and so 
it is not clear how close they have come to M3 1 in the past. As 
discussed by Pawlowski et al. (2011), there is also evidence 
that tidal material can survive in a 'bridge' between the inter- 
acting galaxies for an extended period. Nevertheless, it must 
be cautioned that there is no established precedent for tidal 
dwarfs with the longevity implied by the stellar populations 
of the M31 satellites, and the fact that they continue to adhere 
to such a thin structure is even more perplexing under such a 
scenario. 

There is also another striking characteristic of the observed 



plane. As one will note from examination of Fig. 15 (and in- 
deed the left-hand column of plots in Fig. 10), it is oriented 
perfectly edge-on with respect to the Milky Way. Whilst there 
is a noted bias toward detection of satellites positioned along 
planes oriented in this way, it must be remembered that this 
bias arises primarily due to the propensity for detecting satel- 
lites close to the line of sight passing through M3 1 . Many of 
the satellites observed to lie on our plane are located a good 
distance from this line of sight however and well into the low- 
bias portions of the M31 sky. In any case, the random re- 
alizations of §3.4 suff'er from the same biases and yet show 
unequivocally that the observed plane is very unlikely to arise 
by chance. Hence if we are to accept these results, we must 
also accept the plane's orientation. 

Further to this strikingly edge-on orientation, it is also note- 
worthy that the plane is approximately perpendicular to the 
Milky Way disk. This fact can be easily seen if the con- 
stituent satellites are traced out in Galactic coordinates (i.e. 
all lie on approximately the same Galactic longitude). This 
of course raises the question - how does the orientation of the 
Milky Way's polar plane of satellites compare with this plane? 
Noting that the average pole of the 'Vast Polar Structure' de- 
scribed by Pawlowski, Pflamm-Altenburg, & Kroupa (2012B) 
points roughly in the direction of M3 1 , the two planes are 
approximately orthogonal. These precise alignments are dis- 
cussed in more detail in ILC13, but sufl&ce to say here that 
this alignment is particularly interesting and suggests that the 
Milky Way and M3 1 halos should not necessarily be viewed 
as fully isolated structures. It is entirely conceivable that our 
current ignorance as to the coupling between such structures 
may be to blame for our inability to pin down the precise 
mechanism by which such planes arise. 

Finally, in consequence of these findings, most particularly 
with respect to the highly significant, very thin disk of satel- 
lites that has been identified, it is clear that if ACDM is to 
remain the standard model of cosmology, the occurrence of 
such structures has to be explicable within it. The possible 
deeper implications of the satellite anisotropics are discussed 
in Kroupa et al. (2010), Angus et al. (2011), Fouquet et al. 
(2012) and Kroupa (2012) wherein other alternative cosmolo- 
gies are also highlighted. 



5. CONCLUSIONS 

It is clear that whilst the satellites of M3 1 when taken as a 
whole are no more planar than one can expect from a random 
distribution, a subset consisting of roughly half the sample is 
remarkably planar. The presence of this thin disk of satel- 
lites has been conspicuous throughout the analysis contained 
in this paper. The degree of asymmetry determined from the 
satellite distribution is also found to be relatively high. Of par- 
ticular note, the orientation of the asymmetry is very signifi- 
cant, being aligned very strongly in the direction of the Milky 
Way. When this fact is combined with the apparent orthog- 
onality observed between the Milky Way and M3 1 satellite 
distributions and the Milky Way disk, it appears that the two 
halos may in fact be coupled. Regardless, the great plane of 
satellites identified in this study, and its clear degree of sig- 
nificance, is not directly expected from ACDM cosmology. 
This finding provides strong evidence that thin disks of satel- 
lites do indeed exist in galaxy halos, and whether or not the 
standard model can account for such structures remains to be 
seen. 
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